Where is the string limit in QCD?* 
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The energies of glue in the presence of a static quark-antiquark pair are calculated for separations r ranging 
from 0.1 fm to 4 fm and for various quark-antiquark orientations on the lattice. Our simulations use an improved 
gauge-field action on anisotropic space-time lattices. Discretization errors and finite volume effects are studied. 
We find that the spectrum does not exhibit the expected onset of the universal w/r Goldstone excitations of the 
effective QCD string, even for r as large as 4 fm. Our results cast serious doubts on the validity of treating glue in 
terms of a fluctuating string for r below 2 fm. Retardation effects in the T system are also studied by comparing 
level splittings from the Born-Oppenheimer approximation with those directly obtained in simulations. 



1. Introduction 

Accurate knowledge of the properties of the 
stationary states of glue in the presence of the 
simplest of color sources, that of a static quark 
and antiquark separated by some distance r, is 
an important stepping stone on the way to un- 
derstanding confinement, ft is generally believed 
that at large r, the linearly-growing ground-state 
energy of the glue is the manifestation of the con- 
fining flux whose fluctuations can be described in 
terms of an effective string theory. The lowest- 
lying excitations are then the Goldstone modes 
associated with spontaneously-broken transverse 
translational symmetry. Expectations are less 
clear for small r. 

Even the simplest property, the energy spec- 
trum, of the stationary states of glue interacting 
with a static quark-antiquark pair is not accu- 
rately known. The main goal of this work is to 
remedy this. Here, we present, for the first time, 
a comprehensive determination of the low-lying 
spectrum of gluonic excitations in the presence 
of a static quark-antiquark pair. In this initial 
study, the effects of light quark-antiquark pair 
creation are ignored. A few of the energy levels 
for r less than 1 fm have been studied before [jjj. 
Our results for these quantities have significantly 
improved precision, and we have extended the 
range in r to 4 fm. Most of the energy levels 
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presented here have never been studied before. 
Some of our results were previously reported g|. 

The determination of the energies of glue in 
the presence of a static quark-antiquark pair is 
also the first step in the Born-Oppcnhcimer treat- 
ment of conventional and hybrid heavy-quark 
mesons [|| . The validity of the Born-Oppenheimer 
expansion depends, in part, on the smallness of 
retardation effects. In order to quantify such ef- 
fects for the first time, mass splittings for various 
conventional and hybrid heavy-quark mesons ob- 
tained from the leading Born-Oppenheimer ap- 
proximation are compared to those determined 
in simulations in which the heavy quark prop- 
agates according to a spin-independent nonrela- 
tivistic action. 

2. Computation of the glue energies 

We adopt the standard notation from the 
physics of diatomic molecules and use A to de- 
note the magnitude of the eigenvalue of the pro- 
jection J g ■ f of the total angular momentum J g 
of the gluons onto the molecular axis f . The cap- 
ital Greek letters E,II, A, $, ... are used to in- 
dicate states with A = 0,1,2,3,.. ., respectively. 
The combined operations of charge conjugation 
and spatial inversion about the midpoint between 
the quark and the antiquark is also a symmetry 
and its eigenvalue is denoted by rjcp. States with 
rjcp = 1(— 1) are denoted by the subscripts g (u). 
There is an additional label for the E states; E 



Table 1 

Simulation parameters, including the coupling /?, 
input aspect ratio £, lattice size, and the spa- 
tial link smearing parameters £ and n^. The ap- 
proximate lattice spacings a s , calculated assum- 
ing Tq 1 — 410 MeV, are also given. 



Run 


a s (fm) 


(3 




Lattice 


(C,«c) 


A 


0.29 


2.1 


8 


(10^x20) x 


80 (0.07, 12) 












(0.10, 12) 












(0.15, 12) 


B 


0.29 


2.1 


8 


(10^x20) x 


80 (0.10, 12) 


C 


0.29 


2.1 


8 


14 a x 56 


(0.15, 6) 












(0.10, 12) 


D 


0.27 


2.2 


■5 


12 a x48 


(0.10, 4) 












(0.20, 4) 












(0.30, 4) 


E 


0.22 


2.4 


5 


14 3 x 56 


(0.10, 8) 












(0.15, 8) 












(0.25, 8) 












(0.30, 8) 


F 


0.19 


2.5 5 (10 2 x25) x 


60 (0.10, 12) 












(0.18, 12) 












(0.26, 12) 


G 


0.19 


2.6 


3 


10 a x30 


(0.15, 8) 












(0.30, 8) 


H 


0.12 


3.0 


3 


15 3 x45 


(0.15, 24) 












(0.22, 24) 



states which are even (odd) under a reflection in 
a plane containing the molecular axis are denoted 
by a superscript + (— ). Hence, the low- lying lev- 
els are labelled £+, Ej , £+ E~, E ff , II U , A ff , A„, 
and so on. For convenience, we use T to denote 
these labels in general. 

The glue energies Er(r) were extracted from 
Monte Carlo estimates of generalized Wilson 
loops. Recall that the well-known static poten- 
tial -E s + (r) can be obtained from the large-i be- 
haviour exp[— tE^+ (r)] of the Wilson loop for a 
rectangle of spatial length r and temporal extent 
t. In order to determine the lowest energy in the 
r sector, each of the two spatial segments of the 
r x t rectangular Wilson loop must be replaced 
by a sum of spatial paths, all sharing the same 
starting and terminating sites, which transforms 



as r under all symmetry operations. The easi- 
est way to do this is to start with a single path 
V a , such as a staple, and apply the T projection 
operator which is a weighted sum over all sym- 
metry operations; this yields a single gluon oper- 
ator in the T channel. Different gluon operators 
correspond to different starting paths V a . Using 
several (typically 3 to 22) different such operators 
then produces a matrix of Wilson loop correlators 

Monte Carlo estimates of the Wf 1 (r, t) matrices 
were obtained in eight simulations performed on 
a DEC AlphaStation 500/333 using an improved 
gauge-field action^. The couplings 0, input as- 
pect ratios £, and lattice sizes for each simulation 
are listed in Table [l| Our use of anisotropic lat- 
tices in which the temporal lattice spacing a t was 
much smaller than the spatial spacing a s was cru- 
cial for resolving the glue spectrum, particularly 
for large r. The couplings in the action depend 
not only on the QCD coupling (3, but also on two 
others parameters: the mean temporal link u t and 
the mean spatial link u s . Following Ref. Q, we 
set ut — 1 and obtain u s from the spatial pla- 
quette. We use a s /at — £, the input or bare 
anisotropy, in all of our calculations, accepting 
the small radiative corrections to the anisotropy 
as finite lattice spacing errors which vanish in the 
continuum limit. 

To hasten the onset of asymptotic behaviour, 
iteratively-smeared spatial links Q were used in 
the generalized Wilson loops. A single-link proce- 
dure was used in which each spatial link variable 
Uj (x) on the lattice is mapped into itself plus a 
sum of its four neighbouring (spatial) staples mul- 
tiplied by a weighting factor £. The resulting ma- 
trix is then projected back into SU(3). This map- 
ping is then applied recursively times, forming 
new smeared links out of the previously-obtained 
smeared links. The (£, n^) smearing schemes used 
in the simulations are given in Table |l|. Separate 
measurements were taken for each smearing; cross 
correlations were not determined. The temporal 
segments in the Wilson loops were constructed 
from thermally- averaged links, whenever possi- 
ble, to reduce statistical noise. 

Results for several values of the lattice spacing 



were obtained. Our coarsest lattice (a s ~ 0.29 
fm) was used in runs A, B, and C in order to probe 
very large quark-antiquark separations. A large 
aspect ratio (£ = 8) was needed in order to ade- 
quately resolve the correlation functions. In run 
A, only on-axis measurements were made. Run 
B was done to reduce uncertainties for the large 
f/^s > 9 measurements. To verify the restora- 
tion of rotational symmetry and to rule out prob- 
lems associated with the roughening transition, 
off-axis (r,r,r)/v3 measurements were made in 
run C. Agreement of energies obtained using dif- 
ferent quark-antiquark orientations on the lattice 
also helped identify the continuum A value cor- 
responding to each level (there are only discrete 
symmetries on the lattice). Runs D, E, F, and 
G provided results for different lattice spacings 
a s ~ 0.27,0.22,0.19, and 0.19 fm, respectively. 
Only on-axis measurements were made in these 
runs. Runs F and G correspond to the same 
spatial lattice spacing a s , but have very differ- 
ent temporal spacings at . This provided us with a 
measure of the af errors in our results (our action 
has 0(a^, cQ ) errors). Such information is impor- 
tant for carrying out the a s — > extrapolations. 
Run H provided fine-grained (a s ~ 0.12 fm) mea- 
surements to assist the continuum-limit extrap- 
olations of the and Tl u potentials and was 
useful for verifying the negligible size of quantum 
corrections to the lattice anisotropy. To check the 
anisotropy, the potentials were measured for the 
quark-antiquark axis taken along the very fine- 
grained direction and one of the coarser axes was 
used as the direction of evolution for the sys- 
tem. This run was also important for resolving 
a slight discrepancy between our n„ results and 
those obtained using the simple Wilson ac- 
tion at P = 6.0. 

The matrices (r, t) were reduced in the data 
fitting phase to single correlators and 2x2 corre- 
lation matrices using the variational method. The 
lowest-lying glue energies were then extracted 
from these reduced correlators by fitting a sin- 
gle exponential and a sum of two exponentials, 
the expected asymptotic forms, in various ranges 
twin to i max of the source-sink separation. The 
two-exponential fits were used to check for con- 
sistency with the single-exponential fits, and in 



cases of favourable statistics, to extract the first- 
excited state energy in a given channel. 

Three additional runs on small lattices were 
done to verify that finite-volume errors in our 
results were negligible. We confirmed the small- 
ness of the a s /at renormalization for two values of 
the QCD coupling by extracting the ground state 
potential from Wilson loops in different orienta- 
tions. The hadronic scale parameter ro w 0.5 fm 
was used to determine the lattice spacing 

!■ The 

additive ultraviolet-divergent self-energies of the 
static sources were removed by expressing all of 
our results with respect to E+(r"o). Finite-lattice 
spacing errors were removed by extrapolating our 
simulation results for ro[Er(r) — E s +(ro)] to the 
continuum limit. These extrapolations were car- 
ried out by fitting all of our simulation results to 
an ansatz F cont (r) + -Piatt (?") : a ratio of a poly- 
nomial of degree p+ 1 over a polynomial of degree 
p, where p = 1 or 2, was found to work well for 
the continuum limit form F cont (r), and -Fiatt(f) 
was chosen empirically to be a sum of three terms 
1/y/r, l/r, and 1/r 2 . All fits yielded x 2 /dof near 
unity. Continuum A values were easily identified 
in all cases but one: we were unable to distin- 
guish between a 11^ and <fr u interpretation for the 
on-axis E' u level. 

3. Results 

Our continuum-limit extrapolations are shown 
in Fig. |l|. The ground state £+ is the famil- 
iar static-quark potential. A linearly-rising be- 
haviour dominates the E+ potential once r ex- 
ceeds about 0.5 fm and we find no deviations from 
the linear form up to 4 fm. The lowest-lying ex- 
citation is the n„. There is definite evidence of 
a band structure at large r: the E+', U g , and 
A g form the first band above the ±I U ; the E+, 
E~, W u /$ u , and A u form another band. The E~ 
is the highest level at large r. This band struc- 
ture breaks down as r decreases below 2 fm. In 
particular, two levels, the E~ and E~, drop far 
below their large-r partners as r becomes small. 
Note that for r above 0.5 fm, all of the excitations 
shown are stable with respect to glueball decay. 
As r decreases below 0.5 fm, the excited levels 
eventually become unstable as their gaps above 
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Figure 1. Continuum-limit extrapolations (with uncertainties as indicated) for ro[Er(r)—E- £ + (ro)] against 
r/(2ro) for various T. The dashed lines indicate the locations of the mir/r gaps above the E+ curve for 
m = 1, 2, 3, and 4. The dotted curves are the naive Nambu-Goto energies in four-dimensions. Note that 
we cannot rule out a <!>„ interpretation for the curve labelled 11^. 



the ground state £+ exceed the mass of the light- 
est glueball. 

A feature of any low-energy description of a 
fluctuating flux tube is the presence of Goldstone 
excitations associated with the spontaneously- 
broken transverse translational symmetry. These 
transverse modes have energy separations above 
the ground state given by multiples of n/r (for 
fixed ends). The level orderings and approxi- 
mate degeneracies of the gluon energies at large r 
match, without exception, those expected of the 
Goldstone modes. However, the precise mn/r gap 
behaviour is not observed, as shown in Fig. g. The 
energy differences ro [Er (r) — (r)] for r = II U , 
n g , £~, and T,~ are shown in this figure, along 
with their expected Goldstone mode behaviours, 
indicated by the dashed curves. For separations 
less than 2 fm, one sees from Fig. [I] that the 
gluon energies lie well below the Goldstone ener- 
gies and the Goldstone degeneracies are no longer 
observed. The two E~ states are in violent dis- 



agreement with expectations from a fluctuating 
string. Note also that our results clearly disagree 
with the energies of a Nambu-Goto string naively 
(ignoring quantization difficulties) determined in 
four continuous space-time dimensions. 

These results are rather surprising and cast se- 
rious doubts on the validity of treating glue in 
terms of a fluctuating string for quark-antiquark 
separations less than 2 fm. Note that such a 
conclusion does not contradict the fact that the 
E+(r) energy rises linearly for r as small as 0.5 
fm. A linearly-rising term is not necessarily in- 
dicative of a string: for example, the adiabatic 
bag model predicts a linearly-rising ground state 
much before the onset of string-like behaviour, 
even in the spherical approximation ||. For r 
greater than 2 fm, there are some tantalizing sig- 
natures of Goldstone mode formation, yet signif- 
icant disagreements still remain. To what degree 
these discrepancies can be explained in terms of a 
distortion of the Goldstone mode spectrum aris- 
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Figure 2. Plot of the energy differences ro[Er(r) — 
E^+(r)] against r/(2r ) for T = U u , U g , 
and E~ . Results from all of the simulations are 
shown; the different symbols correspond to the 
different T channels. The pairs of solid curves in- 
dicate the continuum-limit extrapolations. The 
dashed lines show the locations of the mix jr ex- 
citation gaps for m = 1, 2, 3, and 4. 



ing from the spatial fixation of the quark and 
antiquark sources (clamping effect) is currently 
under investigation. This clamping effect can be 
confirmed by comparing our results with the ex- 
citation spectrum of a closed flux in a periodic 
box in which this effect is absent. We also plan 
to study the SU (2) gluon excitation spectrum in 
order to ascertain how much of the spectrum is in- 
dependent of the gauge group. Note that it is very 
unlikely that the transverse extents of our lattices 
are causing the discrepancies with the Goldstone 
spectrum[|j. Results for all r up to 4 fm were 
insensitive to changes in the transverse size from 
0.9 to 3 fm. As an added check, we also calcu- 
lated the wave functions for the first five Fourier 
modes of the naive Nambu-Goto string; we found 
that for r ~ 4 fm, the widths of these wave func- 



tions did not appreciably exceed 1 fm. 

4. Hybrid quarkonium 

Another reason for studying the energies of glue 
in the presence of a quark-antiquark pair is the 
likelihood that these energies will provide insight 
into the nature of hybrid mesons. The study 
of hybrid mesons comprised of heavy quarks is 
the natural starting point in the quest for such 
an understanding. A great advantage in study- 
ing heavy hybrid quarkonium is that such sys- 
tems can be studied not only by direct numerical 
simulation, but also using the Born-Oppcnhcimcr 
expansion. In this approach, the hybrid me- 
son is treated analogous to a diatomic molecule: 
the slow heavy quarks correspond to the nuclei 
and the fast gluon field corresponds to the elec- 
trons^. The first step in the Born-Oppenheimer 
treatment is to determine the energy levels of 
the glue (and light quark-antiquark pairs) as a 
function of the heavy quark-antiquark separation, 
treating the heavy quark and antiquark simply 
as spatially- fixed color sources. Each such en- 
ergy level defines an adiabatic potential. The 
quark motion is then restored by solving the non- 
relativistic Schrodinger equation using these po- 
tentials. Conventional quarkonia arise from the 
lowest-lying potential; hybrid quarkonium states 
emerge from the excited potentials. 

Once the gluon energies are determined, the 
Born-Oppenheimer approach yields the entire 
leading-order spectrum very easily, in contrast to 
direct simulations. However, the validity of the 
Born-Oppenheimer approach relies on the small- 
ness of retardation effects. One way of quan- 
tifying retardation effects is to compare mass 
splittings as determined from the leading Born- 
Oppenheimer approximation with those deter- 
mined from simulations. A Monte Carlo study 
with j3 = 3.0, £ = 3 using a leading-order (but 
lattice-spacing corrected) nonrelativistic action 
for the heavy quark was performed. The effec- 
tive masses for the two lowest-lying mesons in 
the exotic 1 ^ channel are shown in Fig. |[ The 
comparison of mass splittings (for several quark- 
spin-restored states relative to the T) are shown 
in Table 0. All results are expressed in terms of 
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Figure 3. Effective mass plots for the two lowest- 
lying mesons in the exotic T\~ + channel. Mass 
estimates from fitting the 2x2 correlator matrix 
elements using their expected asymptotic forms 
are shown. The heavy quark propagates accord- 
ing to a spin-independent nonrelativistic action. 

the inverse hadronic scale r^ 1 . In the NRQCD 
simulations, the bare quark mass was taken to be 
a s Mb — 2.56. The so-called kinetic mass of the 
T was then determined from its low-momentum 
dispersion relation. Half of this mass was used for 
the quark mass in the leading Born-Oppenheimer 
calculation. This ensured that the T kinetic 
masses were identical in both calculations. As- 
suming that the simulation results do not suf- 
fer significantly from lattice artifacts, one sees 
that retardation affects the spin-averaged mass 
splittings by less than 10%, validating the Born- 
Oppenheimer expansion. 

5. Conclusion 

The spectrum of gluon excitations in the pres- 
ence of a static quark- antiquark pair was compre- 
hensively surveyed for separations r ranging from 
0.1 to 4 fm. Our results raised serious doubts on 
the validity of treating glue in terms of a fluctu- 



Table 2 

Comparison of meson mass splittings as deter- 
mined from the leading Born-Oppenheimer ap- 
proximation (LBO) and nonrelativistic simula- 
tions (NRQCD). The splittings are all taken rela- 
tive to the mass of the T and are in terms of r^ 1 . 
The orbital angular momentum L and glue en- 
ergy r assignments in the LBO for each level are 
indicated, along with the radial quantum num- 
ber n. The size of the differences are listed as 
percentages of the simulation results. 

nLr NRQCD LBO Difference 



Xb 

i-+ 

v -+ 

o+- 
*++ 



2SV+ 

2Pn„ 
1S-- 



0.96(1) 
1.30(1) 
3.29(5) 
4.20(7) 
3.51(8) 
3.56(8) 



0.872(5) 
1.224(3) 
3.166(3) 
3.772(2) 
3.166(3) 
3.807(12) 



9(1)% 
6(1)% 
4(2)% 
10(2)% 
10(2)% 
7(2)% 



ating string for r less than 2 fm. For r between 2 
and 4 fm, some tantalizing signatures of Gold- 
stone mode formation were observed, but dis- 
crepancies still remain. We are currently study- 
ing the role of the spatial fixation or clamping 
of the quark and antiquark sources in distorting 
the Goldstone mode spectrum. Future studies of 
the excitation spectrum of the periodically-closed 
flux arc planned. Lastly, retardation effects in 
quarkonium were found to be sufficiently small to 
validate the Born-Oppcnhcimcr expansion. This 
work was supported by the U.S. DOE, Grant No. 
DE-FG03-97ER40546. 
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